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Abstract 

We study the Fisher model describing natural selection in a popula- 
tion with a diploid structure of a genome by differential- geometric meth- 
ods. For the selection dynamics we introduce an affine connection which is 
shown to be the projectively Euclidean and the equiaffine one. The selec- 
tion dynamics is reformulated similar to the motion of an effective particle 
moving along the geodesic lines in an 'effective external field' of a tensor 
type. An exact solution is found to the Fisher equations for the special 
case of fitness matrix associated to the effect of chromosomal imprinting 
of mammals. Biological sense of the differential- geometric constructions 
is discussed. The affine curvature is considered as a direct consequence of 
an allele coupling in the system. This curving of the selection dynamics 
geometry is related to an inhomogenity of the time flow in the course of 
the selection. 
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1 Introduction 

The selection dynamics in biological populations was usually investigated by the 
methods of the dynamic system theory. Apart from simple cases of explicit 
integrability, the basic problem in such an approach is to find attractors of various 
types and to study the stability problem 0. A survey of basic results of the 
several past decades can be found in monographs || |5|. [T| . 

In common practice of selection dynamics, the geometric methods do not at- 
tract too much attention, although the importance is stressed sometimes |H. M. 
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On the other hand, these methods have demonstrated their efficiency in theo- 
retical and mathematical physics. In modern theoretical physics we observe a 
trend to formulate dynamical principles (having structural resemblance to the 
respective biological relations) in terms of differential geometry and Lie algebra. 
Geometric and algebraic methods provide both suitable mathematical construc- 
tions bringing a system to integrability and they facilitate the study of global 
characteristics of the system. From this standpoint the geometric framework may 
have considerable interest for the mathematical models of biological systems. 

Population models of biological societies display a natural hierarchy regarding 
to the degree and the character of coupling between system components. In the 
most general form, the population dynamics should be considered in the context 
of the ecological system dynamics. The population enters the system as a single 
element connected with the rest ones by trophic, compete and other links. The 
system state is characterized by its number, by genetic structure, by age and sex 
distribution and some other quantities. A complete description of the selection 
dynamics (the change of the population genetic structure) implies to consider 
a multi-locus system with all the genes involved in the selection process. To 
simplify the problem we may consider a single- locus model. A classical example 
of such a case is the well known Fisher model describing a coupling of alleles of the 
given locus in the course of the selection. The Darwin system model suggested 
by Eigen J7J is an utmost case of this simplification. It realizes the basic idea of 
the natural selection: a differential survive of convariantly self-reduplicating units 
(according to Timofeev-Resovskii nomenclature ||), is realized in this model. 

In Ref. we develop Hamiltonian form and thereby simplectic geometric 
description for the selection dynamics in the populations with haploid structure 
of a genome (when any gene is represented in a single instance (single allele)). 

In the geometric framework the evolution of the Darwin system with a stable 
organization is represented as a motion with a constant velocity in an Euclidean 
plane space of information variables. This is a result of absence of a coupling 
between separate genotypes (quasispecies, by Eigen) in the course of the selection. 

In the given work the differential- geometrical methods are applied to the 
Fisher model which is one of the upper level in the above hierarchy with respect to 
the Eigen model. The Fisher model describes the natural selection in populations 
with diploid structure of a genome. In this model, the necessary condition for the 
reproduction (and therefore the selection) of the individuals in the population, 
is a pairwise coupling of separate haploid genotypes (gametes) by the zigote 
formation, i.e. the cells or the organisms with the diploid (double) gene setting. 
The gamete coupling in the population is considered by the Fisher model in terms 
of pairwise coupling of different alleles of the single gene (single-locus model) or 
many genes (multi-locus model) being inherited by the individual from its parents. 

The basic aim of the paper is to clarify the question how the allele coupling 
creates the geometry of the space associated with the selection dynamics. We 
introduce an affine connection related to the first derivative of the Fisher equa- 
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tions in the population variables space. The affine connection is turn out to be 
projectively Euclidean one and simultaneously it is an equiaffine connection. The 
above derivative is presented in the form of equations of motion of an effective 
particle moving along the geodesic lines in an 'effective external field' of a tensor 
type. The Fisher equations are integrated for the fitness matrix of a special form 
associated with the phenomenon of chromosomal imprinting and different ional 
methylation of DNA in the course of the gamete maturation of mammals UTO] . In 



this case the Fisher dynamics is reduced to the dynamics of the Darwin system 
by a suitable replace of time variable. An approximate solution is constructed in 
the weak 'external field'. The equiaffmity is shown to result in the conservation of 
a volume in the population space with respect to the geodesic flows. The volume 
is defined using a certain density of the affine space. 

Mention that the curving of the Euclidean flat space of informational variables 
is the effect of the allele coupling in the system. This curving could be also 
interpreted as an effect of the inhomogenity of the time flow in the in the course 
of the selection. 

2 The Fisher model 

Consider a population where the success of reproductivity (and therefore the se- 
lection dynamics) is determined by a single gene (locus) having N alleles. Denote 
by p a a portion of a- allele in the population, a, (3, . . . = 1, . . . , N. Evidently 

TV 

5>« = i- (i) 

a=l 

Choose the quantities p a as population variables. Then the Fisher system is 
written as ||T|: 

N N 

P a =Pa(J2 UJ o c /3P/3- ^^PfiPu)- (2) 

Here p a = dp a (t)/dt, t is the time. The quantity 

JV 

^a(3Pf3 

is the specific rate of reproduction of a— th allele, and 

N 

is an average specific rate of the population growth as a whole. The matrix uj a p 
(= const) in the population biology is identified to the fitness matrix associated 
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with separate alleles. Let us take the N— th allele as a gauge one and introduce 
the variables 

^ = ln^. (3) 

Pn 

Here and below i,j,k,... = 1, . . . , N — 1. The meaning of the variables z l is that 
they reflect an information quantity per a degree of freedom in the population in 
the frame of Shannon formalism (we mean the information about the genetic 
structure of the population). 



Theorem 1 In the variables (Qj the Fisher system is reduced to the /< 



orm: 



z* = a l + ^^^\ ( 4 ) 



Here, 

Cti = OJiN — ^NN, fiik = Uik — OJNk — UiN + UJNN, (r\ 

n = E k exp(z k ). {0) 

Proof. Introducing the auxiliary variables y a = lnp a , we have from ([!]): 
J2a ex v{ya) = 1- Then exp(y N ) = 1 — Sfc ex P(2/fc)- It is easv to verify that 

( exp(^) exp(V) 

exp(2)= i-&exp W ' expW = TTsP 

E fc exp(y fc ) = y^q> e MVN) = ( 6 ) 
Eqs. (|2]) in the variables t/ q take the form: 

y a = J2 U aP ex P(yf3) + "V eX P(^ + Vv)' ( 7 ) 

p w 

Writing Eqs. (0) for a — i and for a = N and subtracting one from another, we 
obtain (^) using djj). 

3 Projectively Euclidean space 

Let us differentiate the system (|4j) with respect to t and write down the result as 
follows: 

' £i = mO\ E[^ ex P(^) + Suexp{z k )}z k z l + 

1 + (8) 

Introduce the affine connection space .Ajv-i(r^). In the coordinates z 1 the con- 
nection is: 

ni = 2(i + ^) fe ex p(^) + 5 « ex p(^)i- ( 9 ) 
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Eq. (pi) has the form of equations of motion for an effective particle moving along 
the geodesic lines in the affine connection space ^4jv-i(r^) with the 'external 
field' defined by the tensor field A\: 



kl k 



Ai / a x^ln(l + n) 

Al = ( ai + (3 ik ) — ^gpr^-- ( n ) 



where, 



The original equations play the role of constraints to Eqs. (|T 
Consider the basic geometric properties of the space An-i(XIi)- 
The curvature tensor 

f)V q F)T q 
lkL ~ dz k dz l ^ kp H ^ p ki 

can be written using @ as follows: 



«^ = (i + !i ) 1/2 [a^ 7 (i + !i )- 1/2 ^- 
^0 + «r^]. 



The Ricci tensor, R qi = J2k-^kqi ■> * s 



(13) 



R ql = (1+ (1)^(2 -N)^§—. (14) 
Then, we can verify that 

V = 2TjvW-V^ ( 15 ) 

The case N = 2 apparently requires special consideration. If iV > 2, the curvature 
tensor satisfies the criterion of the projectively Euclidean space |L2|, p. 540]. Hence, 
we get 

Theorem 2 The space An-i(XIi) supplied with the affine connection is a 
projectively Euclidean space. 

The form of the affine connection Y l kl (|J) leads to the statement that An-i(^\i) 
is an equiaffine space |13|, §41]. In such a space, a volume exists which is invariant 
under the parallel transport of vectors. The volume element spanned on vectors 
£l k has the form 

V = X/ e n...ijv_i^l 1 • • • Cjv-i 

«i,...,ijv_i 
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which conserves under parallel transport of Ck- Here ei x ,,,i N1 
is completely antisymmetric symbol, ei„.jy-i = 1- In the space Atf-\{T)^}, o 
serves as a fundamental density and is defined by the following condition 

v = ding 

For r kl of the form (§), we obtain: £ s T s ks = <91n(l + 9) N / 2 /dz k . Then, we have 
the theorem true. 

Theorem 3 The fundamental density of the space An-i(^m) with the connection 
§) has the form a = (1 + fi)^ 2 . 

The considered geometry is of a special interest when the geodesical equa- 
tions are exactly integrated. In this case the allele coupling in our approach 
is represented as a pure curving of the population variable space. It is con- 
nected with the following form of the fitness matrix uo a p in the Fisher system (0): 
ujij = oJNjij 7^ N; ^ ujjN, i j. As a variant of a real genetic system result- 
ing in such form of the matrix u a p we can suggest a system with the chromosomal 
imprinting found in the course of the gametogenesis of mammals. According to 
10f1 , the result of the imprinting is that the same allele is included into a genome 



of a zigote in active or in non-active state. The allele activity depends on its 
affiliation to the DNA of male or female gamete. Molecular mechanism of the 
chromosomal imprinting is based on the differential methylation of DNA in the 
course of the gamete maturation. 

For weak 'external field' A\, that takes place under the condition 

oti + Pik = £7ik, e < 1, (16) 

(cKj, f3ik are of the form @) the analysis of the motion can be carried out in two 
steps. First, we consider the motion along the geodesic lines, then we account for 
the 'external field' A\ perturbatively. 



4 Integration of the geodesic equations 

The system flTUp has the form of equations of geodesies of the space ^4at_i(F^) 
under the condition A\ — 0. Eqs. (0), with account of (|||), can be written in 
the form: 

^E^^'=0, (17) 

or 

F + i^ln(l + fi) = 0. 
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The first integration gives 

» - r+sr < 18 > 

where Cj— are integration constants. 

Notice that the original equations (^) under the condition (|T^) take the form: 

a _ a J i_ gSfc7jfcexp(z fc ) 

1 + + ' 1 yj 

Putting e = 0, we see that the restrictions being imposed by Eqs. (f|) on the 
system (|17D, are reduced to 

^ = a 3 . (20) 

Let us now integrate Eqs. fll8|) under the condition ( p0|) . Without loss of gener- 
ality we can assume 

z j = aj <p(t) + b r (21) 

Here bj are the integration constants, and the function (p(t) is determined by the 
equation 

<f>(t) = T^Pq' fi = E^P^W + fofe )' 
whose implicit solution is 

p(t) + ^ «fc 1 exp(a kV (t) +b k )=t- t . (22) 

k 

Let us redefine the parameter in the equations of integral lines ( |2"T| ) setting r = 
(p(t). Then the time t is explicitly expressed in terms of the parameter r as 

t - t = t + ^"fc 1 exp(« fe r + b k ). (23) 

Eqs. ( |2lD take the form of equations of the motion with the constant velocity aj 
with respect to the 'new time' r for an effective particle: 

z^ajr + bj. (24) 

For a weak 'external field' A\ fllBD, Eqs. (|T0f) can be solved approximately as 
follows. 

Let us put 

z^ajr + bjir) (25) 

and take bj to be a slowly varying function of r. Substituting (^) into ( |T§D and 
taking into account (|20"D, we get: 

= e £7i«exp(av + &«)■ ( 26 ) 
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In the first approximation (putting b s to be constants in the right-hand side of 
(|26|) ) we obtain: 

z 3 = ajT + £^7 jfc a A : 1 exp(a fc + 6 fc (0)), (27) 

k 

where 6fc(0) are the constants of integration. 

5 Conclusion 

Discuss a biological sence of the parameters a iy in (|). If /3y = then the 
system (f|) is reduced to the dynamics of haploid populations. A detail analysis of 
such systems was performed in 0. The parameters a, in this case have a meaning 
of specific rates of reproduction of separate genotypes (the Maltusian parameters 
in ecological terms). Thus, it is the matrix (not the complete matrix uj a p) 
which is effectively responsible for the alleles coupling. It is not evident from the 
original Fisher system (0). The allele coupling is represented as a curvature in 
the geometry of the population variable space. The tensor ([12]) is the measure 
of this curvature and it is completely determined by the value ln(l + Q). Note 
that this value is expressed in the population variables in the following form: 
ln(l + ft) — ln(l/pjv)- hi its turn, ln(l/pjv) is a measure of an information of 
the N— th allele in the population by Shannon. It is clear, the more rarely the 
N— th allele occurs in the population the larger is the curvature of the space 
associated with the selection dynamics. Let us also emphasize that the basic 
geometrical characteristic of this space is completely determined by encounter 
frequency of iV-th allele. Among all the possible variants of the matrices u a p, a 
particular interest has the case e = in ([16]) when the system @ is integrated 
explicitly and its solution is reduced to the uniform rectilinear motion (^3) by the 
redefinition of time variable. This case is similar to the Eigen selection dynamics 
in haploid populations and, as it is shown in ||, admits the Hamiltonian form 
of the dynamic equations. Let us note that in this case the allele coupling is 
effectively reduces to time flow change according to (|23|) . 
The work was supported by RFFR grant 98-02-16195. 
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